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Abstract 


We explore the implications of t-b-r (and b-r ) Yukawa coupling unification condi¬ 
tion on the fundamental parameter space and sparticle spectroscopy in the minimal 
gauge mediated supersymmetry breaking (mGMSB) model. We find that this sce¬ 
nario prefers values of the CP-odd Higgs mass > 1 TeV, with all colored sparticle 
masses above 3 TeV. These predictions will be hard to test at LHC13 but they may 
be testable at HE-LHC 33 TeV or a 100 TeV collider. Both t-b-r and b-r Yukawa 
coupling unifications prefer a relatively light gravitino with mass < 30 eV, which 
makes it a candidate hot dark matter particle. However, it cannot account for more 
than 15% of the observed dark matter density. 
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1 Introduction 


The apparent unification of the Standard Model (SM) gauge couplings at the scale Mgut — 
2 x 10 16 GeV in the presence of low scale supersymmetry(SUSY) [1] is nicely consistent with 
supersymmetric grand unified theories (GUTs), but it does not significantly constrain the 
sparticle spectrum. On the other hand, imposing t-b-r Yukawa coupling unification (YU) 
condition at Mgut [2] can place significant constraints on the supersymmetric spectrum in 
order to fit the top, bottom and tau masses [3, 4, 5]. Most work on t-b-r YU condition has 
been performed in the framework of gravity mediation SUSY breaking [6]. Some well known 
choices for GUT-scale boundary conditions on the soft supersymmetry breaking (SSB) 
terms yield, with t-b-r YU condition, quite severe constraints on the sparticle spectrum [3, 
4, 5, 7, 8], which is further exacerbated after the discovery of the SM-like Higgs boson with 
mass nth ~ 125 — 126 GeV [9, 10]. 

Models with gauge mediated SUSY breaking (GMSB) provide a compelling resolution 
of the SUSY flavor problem, since the SSB terms are generated by the flavor blind gauge 
interactions [11, 12, 13]. In both the minimal [11] and general [13] GMSB scenarios, the 
trilinear SSB A-terms are relatively small at the messenger scale, even if an additional sector 
is added to generate the terms [14]. Because of the small A-terms, accommodating 

the light CP-even Higgs boson mass around 125 GeV requires a stop mass in the multi-TeV 
range. This, in turn, pushes the remaining sparticle mass spectrum to values that are out 
of reach of the 14 TeV LHC [15, 16]. There is hope, however, that the predictions can be 
tested at 33 TeV High-Energy LHC (HE-LHC). This tension between the Higgs mass and 
the sparticle spectrum can be relaxed if we assume the existence of low scale vector-like 
particles that provide significant contribution to the CP-even Higgs boson mass [17, 18]. 
We will not consider this possibility here. 

In this paper we investigate YU in the framework of minimal GMSB (mGMSB) sce¬ 
nario. Here, the decoupled messenger fields generate SSB terms of particular non-universal 
pattern. We know from previous studies in gravity mediation SUSY breaking scenarios 
that non-universal SSB terms are necessary to achieve successful YU [5]. Moreover, be¬ 
low Mqut, the mGMSB model contains messenger fields that modify the evolution of the 
gauge and Yukawa couplings, thereby affecting YU. We consider mGMSB in which the 
messenger fields reside in 5 + 5 of 57/(5), which is the simplest scenario that preserves 
unification of the MSSM gauge couplings. We note that previous studies of mGMSB only 
considered evolution below the messenger scale, which is usually much lower than Mgut, 
and therefore could not address the question of YU. 

The remainder of this paper is organized as follows. In Section 2 we briefly outline 
the mGMSB framework that we use for our analysis. We describe our scanning procedure 
along with the experimental constraints we applied in Section 3. In Section 4 we present 
our results, focusing in particular on the sparticle mass spectrum. The table in this section 
presents some benchmark points which summarize the prospects for testing these predic¬ 
tions at the LHC. Our conclusions are presented in Section 5. In Appendix A, we briefly 
discuss the renormalization group equations (RGEs) from the messenger scale to Mgut- 
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2 Minimal Gauge Mediation 


In GMSB models supersymmetry breaking is communicated from a hidden sector to the 
superpartners of SM particles via some messenger fields [12]. In the minimal version, the 
messengers interact with the visible sector only via SM gauge interactions and to the hidden 
sector through an arbitrary Yukawa coupling. The theory is described by the superpotential 

W = Wmssm + (1) 

where Wmssm is the usual superpotential of MSSM, S is the hidden sector gauge singlet 
chiral superfield, and T, <f> are the messenger fields. In order to preserve perturbative gauge 
coupling unification, the messenger fields are usually taken to form complete vector-like 
multiplets of SU(5). Thus one could have N 5 (1 to 4) pairs of 5 + 5, or a single pair of 
10 + 10, or the combination 5 + 5 + 10 + 10. For simplicity, we only consider the case with 
N 5 (one to four) pairs of 5 + 5 vector-like multiplets. 

The singlet field S develops non-zero VEVs for both its scalar and auxiliary components, 
S ^ = (S)+6 2 (F), from the hidden sector dynamics. This results in masses for the bosonic 
and fermionic components of the messenger superfields, 


nib = M n 


1 ± 


A 


Mr, 


rrif = M n 


( 2 ) 


where M mess = A (S) is the messenger scale, and A = (F) / ( S ) sets the scale of the SSB 
terms. Note that (F) < A ( S } 2 to avoid tachyonic messengers, and in many realistic cases 
(F)<\(S) 2 . 

Below the messenger scale M mess , the fields $ and $ decouple generating SSB masses 
for MSSM fields. The masses of the MSSM gauginos are generated at M mess from 1-loop 
diagrams with messenger fields. In the approximation (F) <C A (A) 2 , the gaugino masses 
are given by 

CM • 

Mi(M mess ) = -AiV 5 A, (3) 

471 " 

where i = 1,2,3 stand for U(1)y, SU( 2) l and SU( 3) c sectors, respectively. Since the 
MSSM scalars do not couple directly to the messengers, their SSB masses are induced at 
two loop level: 

mJ(M mess ) = 2N 5 A 2 ^ Q(</>) (^) 2 , (4) 

2=1 

where Ci((j>) is the appropriate quadratic Casimir associated with the MSSM scalar field 


The trilinear A-terms only arise at two-loop level, and thus are vanishing at the messen¬ 
ger scale. However, they are generated below the messenger scale from the RGE evolution. 
The bilinear SSB term also vanishes at M mess , although it is often ignored. We do not 
impose the relation B = 0, anticipating that the value needed to achieve the electroweak 
symmetry breaking (EWSB) can be explained by some other mechanism operating at the 
messenger scale [14]. 

The mGMSB spectrum is therefore completely specified by the following set of param¬ 
eters: 

Mme SS , A, N 5 , tan/3, sign(n) i Cgrav •> (5) 
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where tan/5 is the ratio of the two MSSM Higgs doublet VEVs. The magnitude of /j is 
determined by EWSB condition at the weak scale. The parameter c grav (> 1) affects the 
decay rate of the sparticles into gravitino, but it does not affect the MSSM spectrum, and 
we set it equal to unity. 

It is worth noting that the gaugino masses in mGMSB are non-universal at the mes¬ 
senger scale (if M mess < Mqut)- However, identical relations at the same scale among the 
gauginos can be obtained through RGE evolution in the gravity mediated SUSY breaking 
framework starting with universal gaugino masses at Mqut- On the other hand, Eq. (4) 
implies a universal SSB mass terms for the up (H u ) and down (Ha) Higgs doublets at the 
messenger scale. Within the gravity mediated framework this relation can be obtained 
at the same scale if one imposes non-universal Higgs masses m 2 Hu < rn 2 Hd at Mqut- In¬ 
deed, this non-universality is a necessary condition for compatibility with t-b-r YU in the 
gravity mediation case with universal gaugino masses at Mqut [7, 8]. However, the rela¬ 
tion rn 2 H/ < rm? Hd is introduced by hand in gravity mediated SO(IO) GUT, while t-b-r YU 
in mGMSB can be achieved without such ad hoc conditions, as we will see later. From 
Eq. (4) the sfermion masses in mGMSB are non-universal at the messenger scale, but it is 
not easy to find the same pattern at M mess within the gravity mediation scenario. Hence, 
the sfermion sector could be a good place to see differences between gravity and gauge 
mediation SUSY breaking scenarios that are compatible with YU condition at the GUT 
scale. 


3 Scanning Procedure and Experimental Constraints 

For our scan over the fundamental parameter space of mGMSB, we employed ISAJET 7.84 
package [19] that we modified to include the RGE evolution above M mess in MSSM with 
N 5 pairs of 5 + 5 messengers. In this package, the weak-scale values of gauge and third 
generation Yukawa couplings are evolved from Mz to M mess via the MSSM RGEs in the 
DR regularization scheme. Above M mess the messenger fields are present, which changes 
the beta-functions for the SM gauge and Yukawa couplings. The modified RGEs between 
M mess and Mqut are given in Appendix A. Mqut is defined as the scale at which g\ = g 2 
and is approximately equal to 2 x 10 16 GeV, the same as in MSSM. We do not strictly 
enforce the unification condition = g L — g 2 at Mqut, since a few percent deviation 
from unification can be accounted by unknown GUT-scale threshold corrections [20]. For 
simplicity, we do not include the Dirac neutrino Yukawa coupling in the RGEs, whose 
contribution is expected to be small. 

The SSB terms are induced at the messenger scale and we set them according to ex¬ 
pressions (3) and (4). From M mess the SSB parameters, along with the gauge and Yukawa 
couplings, are evolved down to the weak scale Mz- In the evolution of Yukawa cou¬ 
plings the SUSY threshold corrections [21] are taken into account at the common scale 
Msusy = sJnii L rni R , where m f ~ L and rnf R denote the soft masses of the left and right-handed 
top squarks. The entire parameter set is iteratively run between Mz and Mqut using the 
full 2-loop RGEs until a stable solution is obtained. To better account for leading-log cor¬ 
rections, one-loop step-beta functions are adopted for the gauge and Yukawa couplings, and 
the SSB parameters are extracted from RGEs at multiple scales: at the scale of its own 
mass, rrii = rrq(mj), for the unmixed sparticles, and at the common scale Msusy for the 
mixed ones. The RGE-improved 1-loop effective potential is minimized at Msusy, which 
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effectively accounts for the leading 2-loop corrections. Full 1-loop radiative corrections [21] 
are incorporated for all sparticle masses. 

We have performed random scans over the model parameters (5) in the following range: 


10 3 GeV < 
10 5 GeV < 
40 < 
N 5 = l, 


A < 10 7 GeV, 
M mess < 10 16 GeV, 
tan /? < 60, 

H < 0 . 


( 6 ) 


Varying N 5 from 1 to 4, we find that the low scale sparticle spectrum does not change 
significantly [15]. Thus we present results only for N 5 = 1 case. Regarding the MSSM 
parameter /i, its magnitude but not its sign is determined by the radiative electroweak 
symmetry breaking. In our model we set sign(/i) = — 1. A negative //-term together with 
a positive gaugino mass M 2 gives negative contribution to the muon anomalous magnetic 
moment calculation [22] which would disagree with the measured value for muon g — 2 [23] . 
This, however, is not a problem, since the sparticle spectrum is quite heavy in our scenario. 
As we will see, the smuons are often heavier than 1 TeV. Hence, the SUSY contribution 
to muon g — 2 is not significant, which also is the case in t-b-r YU scenarios [24, 5]. On 
the other hand, a negative sign of // with positive signs for all gauginos helps to realize 
YU. The reason for this is that in order to implement YU, we require significant SUSY 
threshold corrections to the bottom quark Yukawa coupling. One of the dominant finite 
corrections is proportional to gM 3 [25] and it helps to realize Yukawa unification if this 
combination has negative sign [5]. Finally, we employ the current central value for the top 
mass, m t = 173.3 GeV. Our results are not too sensitive to one or two sigma variation of 
m t [27]. 

In scanning the parameter space, we employ the Metropolis-Hastings algorithm as 
described in Ref. [28]. The data points collected all satisfy the requirement of radia¬ 
tive electroweak symmetry breaking (REWSB). After collecting the data, we impose the 
mass bounds on the particles [29] and use the IsaTools package [30, 24] to implement the 
various phenomenological constraints. We successively apply mass bounds including the 
Higgs [9, 10] and gluino masses [31], and the constraints from the rare decay processes 
B s —>■ ji + [32], b —>■ sy [33] and B u —>■ ru T [34], The constraints are summarized below 
in Table 1. Notice that we used wider range for the Higgs boson mass, since there is an 
approximate error of around 2 GeV in the estimate of the mass that largely arises from 
theoretical uncertainties in the calculation of the minimum of the scalar potential, and to 
a lesser extent from experimental uncertainties in the values of rn t and a s . 


123 GeV <m h < 127 GeV 
rrig > 1.5 TeV 

0.8 x 10“ 9 < BR (B s -> /U/O < 6.2 x 10“ 9 (2a) 

2.99 x 10 -4 < BR(6 ->• sy) < 3.87 x 10“ 4 (2a) 

0.15 < B R^Wr T)A f SSM < 2.41 (3a) 

Table 1: Phenomenological constraints implemented in our study. 


Following ref. [28] we define the parameters R t i )r 


and Rb T which quantify t-b-r YU and 
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b-r YU respectively: 


Ma x(y t , y b , y T ) Ma x(y b , y T ) 

Ktbr Mm(y t , y b , y T ) ’ Kbr Min (y b ,y T ) l J 

Thus, R is a useful indicator for Yukawa unification with R < 1.1, for instance, corre¬ 
sponding to Yukawa unification to within 10%, while R = 1.0 denotes ‘perfect’ Yukawa 
unification. 


4 Results 



l°g(A/GeV) log(M mess / GeV) 

Figure 1: Plots in R tbT — A, R tbT — M mess , M mess — A and p — M mess planes. All points are 
consistent with REWSB. Green points also satisfy mass bounds and B-physics constraints. 
Red points are a subset of green and they are compatible with t-b-r YU. In the top panels 
regions below the horizontal line are compatible with YU such that R tbr <1.1 


Figure 1 displays the regions in the parameter space that are compatible with t-b-r 
YU with plots in R tbT — A, R tbT — M mess , M mess — A and // — M mess planes. All points are 
consistent with REWSB. Green points satisfy the mass bounds and B-physics constraints. 
Red points are a subset of green ones and they are compatible with t-b-r YU. In the top 
panels the dashed line indicates the region consistent with YU with R tbr <1.1. The R t i , T —A 
plane shows that the parameter A in mGMSB is constrained by 10 5 GeV < A < 10 6 GeV 
for t-b-r YU with A% T <1.1. As seen from the Rtbr ~ Af mess plane, YU imposes a far 
less stringent constraint on M mess , namely 10 8 GeV < M mess < 10 14 GeV. We also present, 
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tan 0 tan 0 


Figure 2: Plots in R t b T — tan /5 and Ri )T — tan/5 planes. The color coding is the same as in 
Figure 1. 



tan 0 tan 0 


Figure 3: Plots in rrih — tan/5 plane for t-b-r YU (left panel) and b-r YU (right panel). 
The color coding is the same as in Figure 1 except the Higgs mass bound is not applied in 
these panels. 


for additional clarity, the results in M mess — A plane. The MSSM parameter // lies in the 
range —8 TeV < fx < — 3 TeV for compatibility with YU. We see that there are plenty of 
solutions consistent with the mass bounds and constraints from rare B-decays. However, 
only the thin red stripe is compatible with YU. 

The strong impact of t-b-r YU on the fundamental parameter space is somewhat relaxed 
if only b-r YU is imposed. Figure 2 shows the differences between t-b-r and b-r YU cases 
with plots in R t b T — tan (3 and Rb T — tan /5 planes. The color coding is the same as in 
Figure 1. One sees that, t-b-r YU only allows a very narrow range for tan/5, namely 
54 < tan j3 < 60, and perfect t-b-r YU consistent with the experimental constraints can 
be realized when tan /5 ~ 56. On the other hand, the Ri, t — tan /5 plane shows that the 
range for tan /5 is slightly enlarged, 50 < tan /5 < 60, with b-r YU, and perfect b-r YU is 
obtained for 54 < tan /5 < 57. 

A similar comparison based on the range for the SM-like Higgs boson mass is given in 
Figure 3 with plots in nth — tan /5 planes for t-b-r (left panel) and b-r (right panel) YU 
cases. The color coding is the same as in Figure 1, except that the Higgs mass bound is 
not applied in these panels. The left panel shows that demanding 10% or better YU, the 
lightest CP-even Higgs boson mass can be rn h ^ 119 — 125 GeV for t-b-r case, while it can 
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Figure 4: Plots in mg — mg , m„ L — nig, mg 
is the same as in Figure 1. 
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Figure 5: Plots in M mess — mg, A — mg, R t i )r — mg and Ri, T — mg planes. The color coding 
is the same as in Figure 1. 


be as light as 117 GeV in the case of b-r YU, as shown in the right panel. The panels of 
Figure 3 highlight how stringent the Higgs boson mass constraint is on the fundamental 
parameter space of mGMSB compatible with YU. The left panel shows that the 125 GeV 
Higgs boson mass and t-b-r YU more or less equivalently constrain the parameter space. 
On the other hand, as seen from the right panel, b-r YU can be realized for tan/3 > 40, 
while the mass bound on the Higgs boson excludes the region with tan f3 < 50. 

Figure 4 displays the mass spectrum of supersymmetric particles in nig —mg , rriy L — nig, 
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Point 1 

Point 2 

Point 3 

Point 4 

Point 5 

A 

0.94 x 10 6 

1.51 x 10 6 

7.6 x 10 5 

4.7 x 10 5 

1.16 x 10 6 

M mess 

1.5 x 10 14 

1.3 x 10 14 

1.2 x 10 8 

8.6 x 10 13 

1.3 x 10 15 

tan f3 

56.1 

59.3 

58.7 

54.8 

52.7 

h 

-4906 

-6727 

-2995 

-4268 

-3984 

At 

-3564 

-5479 

-1933 

-1834 

-4513 

Ab 

-3600 

-5540 

-1946 

-1891 

-4868 

Ar 

-268.1 

-462.9 

-123.2 

-140.4 

-535.1 

m h 

124 

125 

123 

123 

124 

m H 

1148 

1640 

929 

1665 

2378 

m A 

1141 

1629 

923 

1654 

2363 

m H ± 

1153 

1643 

934.2 

1668 

2380 

rri^o 

Vl,2 

1305, 2438 

2098, 3868 

1045, 1971 

654, 1249 

1602, 2960 

m x° 3A 

4902, 4903 

6725, 6725 

3001, 3002 

4190, 4190 

3993, 3995 

m xG 

2440, 4902 

3871, 6724 

1972, 3002 

1251, 4154 

2961, 3996 

m- g 

6113 

9350 

4995 

3372 

7178 

m UL ] R 

9545, 8590 

13434, 12152 

6903, 6469 

8140, 7206 

7818, 7255 

m h.2 

6325, 7784 

9152, 11112 

5569, 6186 

5011, 6515 

5752, 6745 

m t 

dL,R 

9545, 8317 

13434, 11796 

6903, 6412 

8140, 6937 

7819, 7124 

m b, 2 

6280, 7743 

9134, 11053 

5617, 6148 

5054, 6482 

6022, 6709 


5181 

7122 

2727 

4659 

3560 

mp T 

4843 

6625 

2659 

4358 

3372 

m e L ,R 

5186, 3698 

7129, 5044 

2735, 1478 

4663, 3324 

3567, 2378 

m fl,2 

2600, 4845 

3397, 6628 

1193, 2667 

2340, 4357 

1748, 3381 

Rfbr 

1.00 

1.09 

1.09 

1.06 

1.23 

Rbr 

1.00 

1.08 

1.06 

1.02 

1.07 


Table 2: Benchmark points for exemplifying our results. All masses and scales are given 
in GeV units. Point 1 depicts a solution with perfect YU. Point 2 shows a solution with 
125 GeV mass for h consistent with YU. Point 3 displays a solution with essentially the 
lowest stau mass consistent with YU. Similarly, point 4 shows essentially the lightest mass 
for the lightest neutralino. This point also demonstrates the lightest gluino and stop masses 
obtained from our scans. Point 5 displays a solution which is compatible with b-r YU, but 
not with t-b-T YU. 


rrif. 2 — rrif ] and m-± —m^o planes, with the color coding the same as in Figure 1. As seen from 
the top panels, the colored sparticles are quite heavy. Thus t-b-T YU requires mg > 4 TeV. 
Since the SSB trilinear scalar interaction term A t is relatively small [15], the heavier stop is 
expected to be about the same mass (> 4 TeV) as the lighter stop. Similarly, the squarks 
of the first two families and the gluino have masses in the few TeV range. The rriq L — nig 
plane shows that nig > 3 TeV and nig L > 6 TeV for consistency with t-b-T YU as well as 
the experimental constraints. 

The bottom panels of Figure 4 display the mass spectra for the lightest sparticles which 
are found to be stau, chargino and neutralino. According to the m,f 2 — rn^ panel, f\ can 
be as light as 1 TeV, while nif 2 > 2 TeV. Similarly, the m-± — rn^o plane displays the 
masses of the lightest chargino and neutralino with the latter > 600 GeV, and the former 
> 1200 GeV. 

Finally, Figure 5 displays the gravitino mass, mg , in R t i }r — mg and Ri, t — mg planes 
with the color coding as in Figure 1. As is well-known, the gravitino is usually the LSP 
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in the mGMSB framework, and its mass can be as heavy as 10 TeV. A light gravitino is a 
plausible dark matter candidate and it can also manifest itself through missing energy in 
colliders [36]. 

In the standard scenarios, the relic density bound (Q/r ~ 0.11 [35]) is satisfied with 
a gravitino mass ~ 200 eV [36], which makes it a hot dark matter candidate. The hot 
component of dark matter, however, cannot be more than 15% which, in turn, implies that 
the gravitino mass should be less than about 30 eV [36]. This indeed is predicted from 
Figure 5. The panels of Figure 5 represent the impact of YU on the gravitino mass. Both 
t-b-r and b-r YU do not allow a gravitino that is too heavy, preferring instead a gravitino 
lighter than 100 eV or so. The R t b T — mg plane shows that t-b-r YU is mostly realized 
with mg < 30 eV, with essentially perfect YU possible for nearly massless gravitinos. The 
bound on the gravitino mass is not significantly relaxed for b-r YU. However, essentially 
perfect YU is found for mg < 50 eV. In order to have a complete dark matter scenario one 
could invoke axions as cold dark matter. 

In Table 2 we present five benchmark points that exemplify our results. The points are 
chosen to be consistent with the experimental constraints and YU. The masses are given 
in GeV. Point 1 depicts a solution with perfect YU. Point 2 shows the heaviest mass in our 
scan for the lightest CP-even Higgs boson, while point 3 displays the lightest stau mass 
that we found, consistent with YU. Similarly, point 4 represents the lightest neutralino 
mass foundin our scan. This point also shows the lightest gluino and stop masses that we 
obtained. Point 5 displays a solution which is compatible with b-r YU, but not with t-b-r 
YU. 

5 Conclusion 

Models based on the gauge mediated SUSY breaking are highly motivated, as they provide 
an atractive solution to the SUSY flavor problem. We have explored the implications 
of t-b-r and b-r YU condition on the sparticle spectroscopy in the framework of mGMSB 
where the messengers reside in 5+5 multiplets of SU (5) and only interact with visible sector 
via gauge interactions. We find that YU leads to a heavy sparticle spectrum with the lowest 
mass sparticles being the lightest neutralino (mgi > 600 GeV), chargino (m^± > 1200 GeV) 
and stau (m f > 1000 GeV). In addition, YU prefers values of the CP-odd Higgs boson 
mass mA to be greater than 1 TeV and all colored sparticle masses above 3 TeV. Such 
heavy sparticles can escape detection at the 13-14 TeV LHC, but should be accessible at 
HE-LHC 33 TeV or the proposed 100 TeV collider. We find that YU requirement in the 
mGMSB framework leads to a spectrum of gauginos, Higgses and staus, that is very similar 
to the one in gravity mediated SO(10) GUT, although for different reasons. The details 
of the rest of sfermion spectrum are different, but they are at about the same heavy scale. 
However, YU in mGMSB framework does not require any ad hoc assumptions, such as as 
SSB non-universalities that are necessary in gravity mediated GUTs. It is also interesting 
to note that t-b-r and b-r YU favor a gravitino mass less than 30 eV, which makes it a hot 
dark matter candidate. In this case, the gravitino cannot comprise more than 15% of the 
dark matter density, and one may invoke axions in order to have a complete dark matter 
scenario. 
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A RGEs for MSSM and 5 + 5 


Above the messenger scale M mess the MSSM is supplemented with vector-like messenger 
fields and is described by the superpotential in Eq. (1). Since the theory is supersymmetric 
above M mess , there are no SSB terms in the lagrangian. The N 5 messenger pairs 5 + 5 of 
SU(5) decompose under the SM gauge group as 

3t lf 3 ) ’ 


^ H - $ — L 


M 


i) 2 , - ) + l m 


1) 2, —- ) + D m 


3 + 3 


D 


M 


where numbers in parentheses represent transformation properties under SU( 3) c , SU(2) L 
and U (l)y respectively. 

At the 1-loop level the MSSM gauge couplings contain additional term g^N 5 , where 
a = 1, 2, 3 correspond to the gauge groups U( l)y, SU(2)l and SU(3) C , respectively. At the 
2-loop level additional contributions appear from terms where the Dynkin index is summed 
over all chiral multiplets. Hence the gauge coupling RGEs have the form 


dc/a 

dt 


1 

167T 2 


(A* + siX) 


1 

(I 67 r 2 ) 2 




E 2Ar A?(5K 2 


6=1 



( 9 ) 


Here, tf’ 2) 


are one- and two-loop MSSM beta functions that can be found in Ref. [37] , and 



'7/30 

9/10 

16/15\ 


3/10 

7/2 

° 

( 10 ) 

2/15 

0 

17 / 3 / 



is the contribution from a single 5-plet. Note that 5 and 5 give the same contribution to 
RGEs. 

The anomalous dimensions for MSSM fields are modified only at the 2-loop level, since 
messengers only couple to MSSM fields via gauge interactions. The RGEs for the MSSM 
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Yukawa couplings are 


dyt 

dt 

dyb 

dt 

dy T 

dt 

where /S[ 1,2) are 



one- and two-loop MSSM beta functions [37]. 


( 11 ) 

( 12 ) 

(13) 


References 

[1] H. Baer and X. Tata, Weak Scale Supersymmetry: From Superfields to Scattering 
Events, (Cambridge University Press, 2006); S. P. Martin, Adv. Ser. Direct. High 
Energy Phys. 21, 1 (2010) [hep-ph/9709356]; M. Drees, R. Godbole and P. Roy, Theory 
and phenomenology of sparticles: An account of four-dimensional N=1 supersymmetry 
in high energy physics, (World Scientific,2004). 

[2] B. Ananthanarayan, G. Lazarides and Q. Shafi, Phys. Rev. D 44,1613 (1991) and 
Phys. Lett. B 300, 24 (1993)5; Q. Shafi and B. Ananthanarayan, Trieste HEP 
Cosmol. 1991:233-244. 

[3] H. Baer, S. Raza and Q. Shafi, Phys. Lett. B 712, 250 (2012); A. Anandakrishnan, 

B. C. Bryant, S. Raby and A. Wingerter, Phys. Rev. D 88, 075002 (2013). 

[4] I. Gogoladze, R. Khalid and Q. Shafi, Phys. Rev. D 79, 115004 (2009); S. Antusch and 
M. Spinrath, Phys. Rev. D 79, 095004 (2009); I. Gogoladze, Q. Shafi and C. S. Un, 
JHEP 1207, 055 (2012); I. Gogoladze, Q. Shafi and C. S. Un, JHEP 1208, 028 
(2012); M. Badziak, Mod. Phys. Lett. A 27, 1230020 (2012); A. S. Joshipura and 
K. M. Patel, Phys. Rev. D 86, 035019 (2012); A. Anandakrishnan and S. Raby, Phys. 
Rev. Lett. Ill, 211801 (2013); M. Badziak, M. Olechowski and S. Pokorski, JHEP 
10 , 088 (2013); M. A. Ajaib, I. Gogoladze and Q. Shafi, Phys. Rev. D 88, 095019 
(2013); S. Antusch, S. F. King and M. Spinrath, Phys. Rev. D 89, 055027 (2014); 
M. A. Ajaib, I. Gogoladze, Q. Shafi and C. S. Un, JHEP 1405, 079 (2014); T. Li, 
D. V. Nanopoulos, S. Raza and X. C. Wang, JHEP 1408, 128 (2014). 

[5] I. Gogoladze, R. Khalid, S. Raza and Q. Shafi, JHEP 1012, 055 (2010); M. Adeel 
Ajaib, I. Gogoladze, Q. Shafi and C. S. Un, JHEP 1307, 139 (2013). 

[6] A. Chamseddine, R. Arnowitt and P. Nath, Phys. Rev. Lett. 49 (1982) 970; R. Bar- 
bieri, S. Ferrara and C. Savoy, Phys. Lett. B119 (1982) 343; N. Ohta, Prog. Theor. 
Phys. 70 (1983) 542; L. J. Hall, J. D. Lykken and S. Weinberg, Phys. Rev. D27 (1983) 
2359. 

[7] T. Blazek, R. Dermisek and S. Raby, Phys. Rev. Lett. 88, 111804 (2002). 

[8] H. Baer and J. Ferrandis, Phys. Rev. Lett. 87, 211803 (2001); D. Auto, H. Baer, 

C. Balazs, A. Belyaev, J. Ferrandis and X. Tata, JHEP 0306, 023 (2003). 


11 








[9] G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716, 1 (2012). 

[10] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012). 

[11] M. Dine and A. E. Nelson, Phys. Rev. D 48, 1277 (1993); M. Dine, A. E. Nelson 
and Y. Shirman, Phys. Rev. D 51, 1362 (1995); M. Dine, A. E. Nelson, Y. Nir and 

Y. Shirman, Phys. Rev. D 53, 2658 (1996). 

[12] For a review, see G. F. Giudice and R. Rattazzi, Phys. Rept. 322, 419 (1999). 

[13] P. Meade, N. Seiberg and D. Shih, Prog. Theor. Phys. Suppl. 177, 143 (2009); 
M. Buican, P. Meade, N. Seiberg and D. Shih, JHEP 0903, 016 (2009). 

[14] See, for instance, Z. Komargodski and N. Seiberg, JHEP 0903, 072 (2009) Z. Kang, 
T. Li, T. Liu and J. M. Yang, JHEP 1204, 016 (2012). 

[15] P. Draper, P. Meade, M. Reece and D. Shih, Phys. Rev. D 85, 095007 (2012); 
M. A. Ajaib, I. Gogoladze, F. Nasir and Q. Shah, Phys. Lett. B 713, 462 (2012). 

[16] H. Baer, V. Barger and A. Mustafayev, Phys. Rev. D 85, 075010 (2012) and JHEP 
1205, 091 (2012). 

[17] K. S. Babu, I. Gogoladze, M. U. Rehman and Q. Shah, Phys. Rev. D 78, 055017 
(2008); S. P. Martin, Phys. Rev. D 81, 035004 (2010); P. W. Graham, A. Ismail, 
S. Rajendran and P. Saraswat, Phys. Rev. D 81, 055016 (2010). 

[18] S. P. Martin and J. D. Wells, Phys. Rev. D 86, 035017 (2012). 

[19] ISAJET program, F. E. Paige, S. D. Protopopescu, H. Baer and X. Tata, hep- 
ph/0312045. 

[20] J. Hisano, H. Murayama , and T. Yanagida, Nucl. Phys. B402, 46 (1993); Y. Yamada, 

Z. Phys. C60 (1993) 83; J. L. Chkareuli and I. G. Gogoladze, Phys. Rev. D 58, 055011 
(1998). 

[21] D. M. Pierce, J. A. Bagger, K. T. Matchev, and R.-j. Zhang, Nucl. Phys. B491, 3 
(1997). 

[22] T. Moroi, Phys. Rev. D 53, 6565 (1996) [Erratum-ibid. D 56, 4424 (1997)]. 

[23] G. W. Bennett et al. [Muon (g-2) Collaboration], Phys. Rev. D 80, 052008 (2009). 

[24] H. Baer, C. Balazs, J. Ferrandis and X. Tata, Phys. Rev. D 64, 035004 (2001). 

[25] L. J. Hall, R. Rattazzi and U. Sarid, Phys. Rev. D 50, 7048 (1994); R. Hempfling, Phys. 
Rev. D 49, 6168 (1994); M. S. Carena, M. Olechowski, S. Pokorski and C. E. M. Wag¬ 
ner, Nucl. Phys. B 426, 269 (1994). 

[26] I. Gogoladze, F. Nasir, Q. Shah and C. S. Un, Phys. Rev. D 90, 035008 (2014). 

[27] I. Gogoladze, R. Khalid, S. Raza and Q. Shah, JHEP 1106, 117 (2011). 

[28] G. Belanger, F. Boudjema, A. Pukhov and R. K. Singh, JHEP 0911, 026 (2009); 
H. Baer, S. Kraml, S. Sekmen and H. Summy, JHEP 0803, 056 (2008). 


12 



[29] K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38, 090001 
(2014). 

[30] H. Baer, M. Brhlik, D. Castano and X. Tata, Phys. Rev. D 58, 015007 (1998); 
J. K. Mizukoshi, X. Tata and Y. Wang, Phys. Rev. D 66, 115003 (2002). 

[31] G. Aad et al. [ATLAS Collaboration], JHEP 1410, 24 (2014), ATLAS-CONF-2013- 
062 and JHEP 1409, 103 (2014); S. Chatrchyan et al. [CMS Collaboration], JHEP 
1210, 018 (2012) and CMS-PAS-SUS-13-019. 

[32] R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 110, 021801 (2013). 

[33] Y. Amhis et al. [Heavy Flavor Averaging Group Collaboration], arXiv: 1207.1158 [hep- 
ex], 

[34] D. Asner et al. [Heavy Flavor Averaging Group Collaboration], arXiv: 1010.1589 [hep- 
ex], 

[35] G. Hinshaw et al. [WMAP Collaboration], Astrophys. J. Suppl. 208, 19 (2013). 

[36] M. Viel, J. Lesgourgues, M. G. Haehnelt, S. Matarrese and A. Riotto, Phys. Rev. D 
71, 063534 (2005); J. L. Feng, M. Kamionkowski and S. K. Lee, Phys. Rev. D 82, 
015012 (2010). 

[37] S. P. Martin and M. T. Vaughn, Phys. Lett. B 318, 331 (1993). 


13 



